/ 

AD-A114  160  OPTIMCTRICS  INC  ANN  ARBOR  Ml  ft%  6/1 

LASER  PROPAGATION  RESEARCH.  VOLUME  I.  SLANT  PATH  TRANSMISSION  M— ETC (U) 
NOV  60  W  M  tUTMANr  $  T  HANLEY  OAAD07-SO-C-0056 


PHOTOGRAPH  THIS  SHEET 


LEVEL 


m 

INVENTORY 


Laser  "tVcpacfoLor)  i^esecvzh.  Volume  J 

DOCUMENT  IDENTIFICATION  Mp.  OMI-&0  -OCfc  (jf) 

Q?^Wcf  DftfWO?  -%0-  Q.-QOSC> _ a/qv.  BO 


DISTRIBUTION  STATEMENT  A  I 

Approved  lor  public  release;  | 
Distribution  Unlimited  ! 

DISTRIBUTION  STATEMENT 


ACCESSION  FOR 


NT  IS  GRAAI 

OTIC  TAB 

UNANNOUNCED 


BY 


DISTRIBUTION  / 


AVAILABILITY  CODES 


DIST  AVAIL  AND/OR  SPECIAL 


DATE  ACCESSIONED 


copf 

.  INSPECTED y 


DISTRIBUTION  STAMP 


DATE  RECEIVED  IN  DTIC 

PHOTOGRAPH  THIS  SHEET  AND  RETURN  TO  DTIC-DDA-2 


DOCUMENT  PROCESSING  SHEET 


adA114i6q 


LASER  PROPAGATION  RESEARCH 


OMI-8O-OO8 (I) 


VOLUME  I 

SLANT  PATH  TRANSMISSION  MEASUREMENTS 


William  M.  Gutman 
Stephen  T.  Hanley 


Prepared  for 

A-tmospheric  Sciences  Laboratory 
White  Sands  Missile  Range,  New  Mexico 


Contract  DAAD07-80-C-0056 


November  1980 


OpriMETRi'cs,  Inc. 

2000  Hogback  Rd.,  Suite  3 
Ann  Arbor,  Ml  48104 


distribution  statement  a 

Approved  for  public  release; 
Distribution  Unlimited  i 


82 


05  03  042 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Intend) 


REPORT  DOCUMENTATION  PAGE 


I  REPORT  NUMBER 


7  AUTHOR  (s) 

William  M.  Gutman 
Stephan  T.  Hanley 


9  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

OptiMetrics,  Inc. 

2000  Hogback  Road,  Suite  3 


II  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Atmospheric  Sciences  Laboratory 
US  Army  Electronics  R&D  Command 
WSMR,  NM  88002 


14  MONITORING  AGENCY  NAME  AND  ADDRESS 
(if  different  front  Controlling  Off  at J 


16  DISTRIBUTION  STATEMENT^/  thn  Report 


READ  INSTRUCTIONS 
HFKMU  COMPUTING  FORM 


RECIPIENT  S  CATAL 


5  TYPE  OF  REPORT  ft  PERIOD  COVERED 

Technical  Report 
3/80  -  10/80 


6  PERFORMING  ORG  REPORT  NUMBER 

0M1 -80-008 (I) 


8  CONTRACT  OR  GRANT  NUMBERjj^ 

DAAD07-80-C-0056 


10  PROGRAM  ELEMENT  PROJECT  TASK 
AREA  6  WORK  UNIT  NUMBERS 


Volume  1 
Task  F.3.2 


13  NUMBER  OF  PAGES 

50 


15  SECURITY  CLASS  (of  thn  report) 

Unclassified 

Ts 4  declassification/ DOWNGRADING 
SCHEDULE 


distribution  statement  a 


Approved  lot  public  ieleaBe) 
Distribution  Unlimited 


t7  Distribution  STATEMENT  ^ of  the  abstract  entered  m  Block  20,  if  different  from  Report) 


19  KEY  WORDS  f Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 

High  Energy  Laser  H^O,  HD0,  N^O,  0^,  THC,  CH^ 

Atmospheric  Transmission  _  .  _  r 

Fourier  Transform  Spectrometer  (FTS) 

ARKY  Site 
NOP  Site 


Solar  Spectrum 


20  ABSTRACT  0n  reverse  side  tf  necessary  and  identify  by  block  number) 

An  existing  Fourier  Transform  Spectrometer  (FTS)  data  base  of  slant  path 
atmospheric  transmission  spectra  using  the  sun  as  a  source  is  described.  Re¬ 
cent  improvements  to  the  FTS  measurements  system  are  discussed  and  results  of 
pyroelectric  radiometer  measurements  of  slant  path  solar  transmission  are 
reviewed . 

Slant  path  atmospheric  transmission  modeling  developments  and  use  of  these 
odels  in  the  interpretation  of  slant  path  FTS  transmission  data  are  discussed. 


1  JAN  73  M73  EDITION  OF  ,  NOV  65  IS  OBSOLETE  11  N Q T  .  A  R  R  T  F  T  F.  D 

SECURITY  CLASSIFICATION  OF  THIS  PAGEfIT'ior  Data  t  It  ter  ej) 


UNCLASSIFIED 


security  classification  OF  This  PAGE  (V  hr*  DalaLmlrrrd) _ 

Techniques  for  atmospheric  absorbing  gas  vertical  density  profile  extrac¬ 
tion  are  presented  and  results  of  model  comparisons  with  FTS  data  are 
reviewed. 

Recent  measurements  of  l^O,  N-0,  THC,  CH-,  0^,  and  CC>2  atmospheric 
concentrations  at  the  ARKY  and  NOP  sites  at  WSMR  are  presented  and 
discussed.  Results  of  aerosol  particulate  distribution  and  mass  loading 
measurements  at  the  WSMR-NOP  sites  are  presented  and  described. 


SSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (W'bftl  Dal.i  Intrrtd) 


CONTENTS 


1.  INTRODUCTION .  1 

2.  SOLAR  SPECTRAL  MEASUREMENTS .  2 

2.1  Augmentation  of  Existing  Fourier 

Transform  Spectrometer  Data  Base .  2 

2.2  Investigation  of  the  Use  of  Pyroelectric 
Radiometers  for  Solar  Spectral  Measurements....  2 

2.3  Additional  Instrumental  Capabilities .  5 

2.3.1  Implementation  of  Improved 

FTS  Software .  5 

2.3.2  Installation  of  the  FTS  System  in 

the  Field  Measurement  Van . 8 

3.  ANALYSIS  AND  INTERPRETATION  OF  PREVIOUSLY 

COLLECTED  FOURIER  TRANSFORM  SPECTROMETER  DATA .  10 

3.1  Correlation  of  CO2  and  DF  Laser  Line 

Extinction  Coefficients  with  Atmospheric 
Absorption  Line  Measurements .  10 

3.2  Slant  Path  Transmission  Model  Development .  13 

3.3  Analysis  of  Measured  FTS  Spectra .  22 

3.3.1  Comparisons  of  Slant  Path 
Transmission  Predictions  to 

Measured  FTS  Data .  22 

3.3.2  Extraction  of  Vertical  Density 
Profiles  of  Absorbing  Species 

from  Measured  FTS  Data .  30 

3.3.3  Further  Utilization  of  Existing 

Slant  Path  FTS  Solar  Spectra .  44 

4.  CONCLUSIONS  AND  RECOMMENDATIONS .  4  6 

REFERENCES .  4  9 


11 


LIST  OF  ILLUSTRATIONS 


1.  Relative  Energy  Received  by  the  Detector  During  a 

Sampling  Interval .  4 

2.  Extinction  Coefficient  at  DF  Laser  Line  P,(12) 

at  2445.35  cm-1 .  11 

3.  Extinction  Coefficient  at  DF  Laser  Line  Po(7) 

at  2570.52  cm_l .  12 

4.  Extinction  Coefficient  at  CO2  Laser  Line  P(56) 

at  907.78  cm  .  14 

5.  Extinction  Coefficient  at  CO2  Laser  Line  P(20) 

at  944.19  cm~l .  15 

6.  Extinction  Coefficient  at  C02  Laser  Line  P(12) 

at  9  51.19  cm-l .  16 

7.  Extinction  Coefficient  at  CO2  Laser  Line  R(28) 

at  1083.48  cm”1 .  17 

8.  Extinction  Coefficient  at  the  CO2  Laser  Line 

R  ( 48  )  at  1093.89  cm"! .  18 

9.  Slant  Path  Transmission  Spectra  at  4  ym 

for  87.5°  from  Zenith .  23 

10.  Slant  Path  Transmission  Spectra  at  4  ym 

for  84.4°  from  Zenith .  24 

11.  Slant  Path  Transmission  Spectra  at  4  ym 

for  78.9°  from  Zenith .  25 

12.  Slant  Path  Transmission  Spectra  at  4  ym 

for  78.9°  from  Zenith .  26 

13.  Slant  Path  Transmission  Spectra  at  4  ym 

with  same  Conditions  as  Figure  11 .  27 

14.  Slant  Path  Transmission  Spectra  at  10.6  ym 

for  62.8°  from  Zenith .  28 

15a.  Slant  Path  Transmission .  32 

15b.  Vertical  Profile .  33 

iii 


16a.  Slant  Path  Transmission .  34 

16b.  1^0  Vertical  Profile .  35 

17a.  Slant  Path  Transmission .  36 

17b.  H2O  Vertical  Profile .  37 

18a.  Slant  Path  Transmission .  38 

18b.  H2O  Vertical  Profile .  39 

19a.  Slant  Path  Transmission .  40 

19b.  H2O  Vertical  Profile .  41 

20a.  Slant  Path  Transmission .  42 

20b.  H2O  Vertical  Profile .  43 


1 

INTRODUCTION 

This  report  describes  work  performed  for  the  Atmos¬ 
pheric  Sciences  Laboratory  (ASL)  in  support  of  their  High 
Energy  Laser  Meteorology  (HELMET)  responsibility  to  the 
High  Energy  Laser  Systems  Test  Facility  (HELSTF) .  The 
measurements  and  analyses  reported  herein  continue  and 
extend  work  begun  in  1977.  It  describes  work  completed 
during  the  period  1  March  through  30  September  1980. 

In  the  interest  of  clarity,  this  report  has  been  se¬ 
parated  into  two  volumes.  Volume  1  describes  progress, 
results,  and  data  analyses  associated  with  the  measurements 
at  atmospheric  transmission  spectra  using  the  sun  as  a 
source.  Volume  II  describes  recent  work  carried  out  in 
atmospheric  gas  analysis  studies  and  contains  results  of 
measurements  performed  during  the  most  recent  reporting 
period.  This  volume  is  devoted  to  the  FTS  Solar  Measure¬ 
ments. 
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SOLAR  SPFCTRAL  MEASUREMENTS 


2.1  AUGMENTATION  OF  EXISTING  FOURIER  TRANSFORM 

SPECTROMETER  DATA  BASE 

The  techniques  and  procedures  used  in  the  accumulation 
of  the  data  base  consisting  of  sets  of  solar  source  atmos¬ 
pheric  spectra,  and  the  extraction  of  extinction  coefficients 
from  these  spectra  has  been  reported  previously  [1,2,3]  and 
will  not  be  discussed  in  great  detail.  Measurements  were 
made  on  twelve  days  in  the  first  quarter  of  the  reporting 
period  using  the  original  software  as  supplied  by  the 
Fourier  Transform  Spectrometer  (FTS )  manufacturer  (MIDAC 
Corporation) .  In  the  final  quarter,  trial  measurements  were 
made  on  three  additional  days  using  improved  operating  soft¬ 
ware  supplied  by  MIDAC.  These  measurements  were  made  pri¬ 
marily  as  a  test  of  the  usefulness  of  the  new  system,  but 
the  results  are  available.  Since  a  total  of  only  15  days 
of  measurements  were  made  during  the  present  reporting 
period,  and  since  these  days  are  not  necessarily  representa¬ 
tive  of  each  seasonal  period,  average  extinction  coeffi¬ 
cients  are  not  very  meaningful  except  when  viewed  in  the 
context  of  the  entire  2-1/2  year  data  base.  For  this 
reason,  tabulations  and  plots  of  these  values  are  not  pre¬ 
sented  in  this  report. 

2.2  INVESTIGATION  OF  THE  USE  OF  PYROELECTRIC  RADIOMETERS 

FOR  SOLAR  SPECTRAL  MEASUREMENTS 

During  several  of  the  FTS  data  runs,  two  Laser  Pre¬ 
cision  Corporation  electrically  calibrated  pyroelectric 
radiometers  were  operated  simultaneously  with  the  Fourier 
Transform  Spectrometer.  These  devices  were  provided  by  the 
Air  Force  Weapons  Laboratory,  and  were  operated  in  the  same 


beam  of  solar  radiation  as  the  spectrometer.  Each  was 
equipped  with  a  narrow  bandpass  filter  with  a  transmission 
width  of  about  10  cm  ^  centered  at  about  2585  cm  .  This 
is  within  the  DF  laser  region  of  the  spectrum.  The  majority 
of  the  spectra  collected  with  the  Fourier  Transform  Spectro¬ 
meter  were  collected  at  a  nominal  resolution  of  0.12  cm 
This  value  represents  the  spacing  between  the  first  order 
zero  crossings  of  the  unapodized  sin  x/x  instrumental  line 
shape.  Comparing  a  given  data  point  in  several  spectra, 
the  values  are  proportional  to  the  energy  within  the  resolu¬ 
tion  element  falling  on  the  detector  during  the  scan.  A 
plot  of  the  log  of  this  value  versus  airmass  yields  the 
extinctance  which  is  the  slope  of  this  graph.  The  techni¬ 
que  is  discussed  fully  in  References  1  and  2.  It  had  been 
hoped  that  the  10  cm  1  bandpass  filter  would  be  narrow 
enough  that  a  similar  technique  applied  to  the  radiometer 
data  would  yield  average  extinction  coefficients  over  the 
spectral  region  which  would  correlate  with  those  obtained 
at  high  resolution  with  the  spectrometer.  It  would  then 
have  been  possible  to  monitor  variations  in  extinction 
values  with  much  simpler  and  less  expensive  instrumentation. 

The  pyroelectric  devices  proved  not  to  be  practical 
for  such  use,  however.  They  exhibit  a  long  term  radio- 
metric  equivalent  drift  on  the  order  of  2  microwatts.  The 
solar  flux  transmitted  through  the  bandpass  filter  however, 
is  on  the  order  of  only  a  few  microwatts.  The  radiometer 
response  is  sampled  at  the  rate  of  about  3  Hz.  Figure  1 
is  a  plot  of  the  log  of  the  relative  energy  falling  on  the 
pyroelectric  detectors  during  one  sampling  interval  as  a 
function  of  air  mass  for  each  of  the  two  radiometers  for  a 
typical  day.  Included  for  comparison  is  a  typical  plot 
taken  from  0.12  cm-'*'  resolution  Fourier  transform  spectra 
for  the  DF  laser  line  P-.(6)  frequency  of  2594.20  cm  . 
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FIGURE  1.  RELATIVE  ENERGY  RECEIVED  BY  THE  DETECTOR  DURING  A  SAMPLING  INTERVAL  IS 
PLOTTED  S EM I- LOGARITHMICALLY  VERSUS  AIR  MASS  FOR  THE  TWO  PYROELECTRIC 
RADIOMETERS,  AND  FOR  THE  FTS  SPECTRA  AT  DF  LASER  LINE  P3(6)  FREQUENCY 
AT  2  594.20  cm  THE  SLOPES  OF  THE  CURVES  ARE  PROPORTIONAL  TO  THE 

EXTINCTION  COEFFICIENT. 


The  extent  of  the  error  bars  represents  the*  maximum  and 
minimum  values  observed  dur inc;  a  15-seconu  time  period.  It 
is  clear  that  a  least  squares  fit  to  these  data  would  exhibit 
a  large  standard  deviation,  and  that  the  slopes  are  not  con¬ 
sistent  with  each  other.  In  fact,  in  one  case,  the  energy 
decreases  with  decreasing  air  mass.  It  is  thus  apparent 
that  the  electrically  calibrated  radiometers  are  not  useful 
for  making  measurements  with  such  a  narrow  landpass  filter. 

In  order  to  reduce  the  drift  problem,  the  energy  falling 
on  the  detector  should  be  large  compared  with  the  drift, 
and  this  could  be  achieved  by  using  a  much  wider  bandpass. 

Such  a  measurement,  however,  which  would  involve  a  region 
containing  perhaps  several  dozer,  strong  spectral  lines, 
would  bear  little  relation  with  the  high  resolution  spectral 
measurements . 

2.3  ADDITIONAL  INSTRUMENTAL  CAPABILITIES 

2.3.1  IMPLEMENTATION  OF  IMPROVED  FTS  SOFTWARE 

A  major  part  of  the  FY  1980  effort  was  devoted  to  the 
implementation  of  the  improved  MIDAC  Corporation  software 
for  the  spectrometer  system.  The  primary  motivation  for  the 
acquisition  of  the  software  update  was  to  be  able  to  perform 
Fourier  transforms  on  larger  data  sets.  The  previous  system 
was  limited  to  transforms  of  131072  points,  and  in  order  to 
eliminate  aliasing  problems,  this  constraint  effectively 
limited  the  resolution  to  0.12  cm  ,  where  the  resolution 
criterion  is  the  spacing  between  the  first  order  zero  cros¬ 
sings  of  the  unapodized  sin  x/x  instrument  lineshape.  This 
resolution  proved  to  be  insufficient  to  answer  subtle  ques¬ 
tions  about  the  shapes  of  atmospheric  CO^  absorption  lines. 

MIDAC  Corporation  chose  to  completely  revise  its  spectro¬ 
meter  operating  software.  The  foremost  improvement  was  the 
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doubling  of  the  size  of  the  data  set  which  could  be  trans¬ 
formed.  This  of  course  doubled  the  resolution  which  could 
be  used  without  aliasing  problems.  Another  significant 
improvement  was  the  utilization  of  the  stepped  gain  feature. 
This  feature,  which  had  always  been  available  in  hardware, 
was  made  accessible  from  software.  In  effect,  after  a  spe¬ 
cified  number  of  data  points,  the  overall  inter fk. ^gram 
signal  gain  before  digitization  is  increased  by  a  factor 
of  8.  The  step  is  made  between  data  points,  and  at  a  loca¬ 
tion  in  the  interferogram  sufficiently  past  the  center  burst 
where  the  dynamic  range  of  *the  signal  is  much  smaller.  The 
net  effect  is  improved  utilization  of  the  analoa-to-digital 
converter  and  improved  signal  to  noise  ratio.  Another  im¬ 
portant  improvement  entailed  full  use  of  a  file  header 
containing  all  pertinent  information  about  each  spectrum.. 
Included  in  this  information  is  the  bit  shift  applied  to 
each  interferogram  before  transformation.  This  shift  is 
applied  to  eliminate  some  noise  bits  prior  to  transforma¬ 
tion,  and  in  effect  scales  the  resulting  spectrum  by  powers 
of  two.  Having  this  information  available  allows  the  direct 
comparison  of  spectra  without  having  to  surmise  the  relative 
scaling  as  was  necessary  under  the  old  operating  system. 
Finally,  the  new  MIDAC  software  is  much  simpler  to  use  for 
the  operator.  The  desired  spectrometer  operating  conditions 
are  entered  directly  into  a  control  file  by  the  operator 
rather  than  by  coding  in  the  form  of  control  words.  The 
net  result  is  that  the  system  can,  in  principle,  be  fully 
utilized  by  an  operator  with  little  training  in  minicomputers 
and  the  intricacies  of  the  hardware. 

V 

Of  course  these  improvements  were  not  made  without  some 
sacrifices.  Under  the  previous  system,  the  interf erograms 
were  collected  into  double  16-bit  integers,  and  the  spectra 
resulting  were  in  the  form  of  single  16-bit  integers.  In 
order  to  eliminate  the  resulting  file  size  discrepancy, 
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and  to  increase  the  dynamic  range  available  in  the  spectra, 
the  new  system  uses  double  words  for  the  spectra  as  well 
as  the  inter ferograms .  This  fact  coupled  with  doubling  the 
size  of  the  data  sets  in  order  to  double  the  resolution 
reduces  the  number  of  spectra  which  can  be  stored  on  a  re¬ 
movable  disc  cartridge  to  two,  a  number  insufficient  for 
calculation  of  extinction  coefficients.  In  addition,  the 
time  required  for  the  Fourier  transformation  of  a  262144 
point  data  set  is  much  longer,  due  both  to  its  size  and  to 
the  fact  that  the  operation  is  done  in  double  precision. 
Rather  than  having  the  capability  to  collect  and  transform 
four  spectra  per  hour  as  before,  the  limit  was  about  one 
per  hour.  MID AC  was  able  to  overcome  the  storage  and  time 
limitations  through  the  use  of  the  magnetic  tape  drive.  A 
modified  version  of  the  software  was  provided  which  can  be 
used  to  store  each  complete  interf erogram  on  a  reel  of 
magnetic  tape  as  it  is  collected.  A  10-inch  reel  of  tape 
can  contain  eleven  high  resolution  data  sets.  After  a  reel 
of  tape  is  filled  with  interf erograms ,  each  file  can  be 
read  from  the  tape  by  the  data  system,  Fourier  transformed, 
and  the  spectrum  written  back  onto  the  tape  over  the  inter- 
ferogram.  This  can  be  accomplished  after  all  data  collec¬ 
tion  is  complete,  and  hence  transformation  does  not  place 
time  constraints  on  data  collection. 

The  data  tapes  which  result  from  the  procedure  described 
above  are  of  course  incompatible  with  those  which  were 
created  using  the  old  operating  software,  so  it  was  neces¬ 
sary  to  modify  all  previously  implemented  data  reduction 
programs  used  with  the  HP-1000  system.  The  double  word 
integers  could  not  be  read  directly,  so  it  was  necessary 
to  construct  a  routine  to  read  the  two  words  separately, 
but  properly  recover  the  double  word  value.  It  was  found 
that  the  values  could  be  converted  to  two  word  real  numbers 
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in  the  Hewlett-Packard  format  without  loss  of  precision. 

\ 

Tne  previously  used  tape  read  subroutine  was  modified  in 
this  manner.  The  point  spacing  in  the  spectra,  of  course, 
was  also  different.  For  archival  storage,  a  further  refine¬ 
ment  was  made.  The  MIDAC  software  created  tapes  consist  of 
512  word  o\r  256  data  point  records  on  800  bpi  9-track  tapes. 
By  writing \a  new  tape  consisting  of  8192  word  or  4096  point 
records,  most  of  the  space  required  for  inter-record  gaps 
could  be  saved.  Furthermore,  it  is  possible  to  save  only 
those  spectral  regions  of  interest,  and  since  the  new  tape 
is  written  on  a  1600  bpi  tape  drive,  the  net  result  is  that 
several  days  of  spectra  can  be  stored  on  a  given  cape.  New 
versions  of  the  data  reduction  tapes  were  then  written  so 
that  these  archival  tapes  can  be  used  directly  in  the  pro¬ 
cedure  . 

Some  difficulties  were  encountered  with  the  new  MIDAC 
software,  initially.  The  ASL  spectrometer  system  utilizes 
a  Versatec  plotter/printer  for  spectral  plotting  which  is 
different  from  that  in  the  MIDAC  test  system.  Several  iter¬ 
ations  were  required  in  the  plotting  routine  to  correctly 
interface  with  the  device,  and  to  produce  plots  with  a  cor¬ 
rect  frequency  scale.  After  the  problems  were  overcome, 
the  system  was  tested  by  obtaining  spectra  on  three  days, 
and  calculating  extinction  values  in  the  usual  manner.  When 
a  workable  routine  had  been  devised,  and  the  system  was 
operating  properly,  the  decision  was  made  to  move  the  spec¬ 
trometer  into  its  van. 

2.3.2  INSTALLATION  OF  THE  FTS  SYSTEM  IN  THE  FIELD 
MEASUREMENT  VAN 

The  van  to  be  used  for  field  measurements  had  been 
previously  prepared  prior  to  the  FTS  system  installation. 
Every  effort  was  made  to  insure  a  weather  tight  and  dust 
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free  environment  for  the  system.  The  roof  was  repaired,  the 
floor  tiled,  the  walls  lined,  and  the  exterior  painted  wTith 
a  highly  reflective  white  paint  to  minimize  the  summertime 
thermal  load.  However,  due  to  uncontrollable  circumstances 
during  the  move,  the  germanium  on  potassium  bromide  beam¬ 
splitter  which  had  been  used  during  the  collection  of  all 
previous  solar  spectra  was  damaged.  A  replacement  was 
ordered,  and  was  being  installed  and  tested  as  the  reporting 
period  ended.  In  the  van,  the  system  will  be  used  to  col¬ 
lect  slant  path  solar  spectra  in  much  the  same  way  as  has 
been  done  previously.  A  Carson  Astronomical  Instruments 
6-inch  heliostat  is  available  for  solar  tracking.  How¬ 
ever,  this  device  required  minor  modification  for  proper 
operation.  The  tracking  system  contains  a  drive  motor  for 
the  right  ascension  axis  and  slew  motors  for  both  the  right 
ascension  and  declination  axes.  When  the  right  ascension 
drive  motor  was  in  operation,  the  slew  motor,  which  should 
have  remained  fixed,  tended  to  rotate  with  the  net  result 
that  the  tracking  mirror  remained  fixed.  After  numerous 
attempts  to  correct  this  problem  by  balancing  the  frictional 
forces  with  a  crude  spring  clutch  incorporated  into  the 
system,  it  became  apparent  that  when  the  frictional  forces 
were  large  enough  to  prevent  the  slew  motor  from  freely 
rotating,  then  the  drive  motor  torque  was  insufficient  to 
overcome  the  same  frictional  forces.  The  situation  was 
remedied  by  modifying  the  electrical  circuit  so  that  a  DC 
voltage  is  applied  across  the  slew  motor  whenever  it  is 
not  being  used  to  slew.  A  full  wave  rectified  85  volt  rms 
AC  voltage  proved  sufficient  to  solve  the  problem. 
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ANALYSIS  AND  INTERPRETATION  OF  PREVIOUSLY  COLLECTED 
FOURIER  TRANSFORM  SPECTROMETER  DATA 

Analysis  of  previously  collected  data  carried  out  during 
this  reporting  period  can  be  divided  into  two  categories: 
a)  attempts  to  correlate  extinction  coefficient  variation 
with  variation  in  other  measurable  parameters;  and  b) 
attempts  to  use  accepted  atmospheric  models  to  predict 
measured  slant  path  transmission.  Additionally,  preliminary 
efforts  to  extract  vertical  profiles  of  infrared  absorbing 
species  based  solely  on  measured  slant  path  transmission 
spectra  have  proved  successful. 

3.1  CORRELATION  OF  C02  AND  DF  LASER  LINE  EXTINCTION 
COEFFICIENTS  WITH  ATMOSPHERIC  ABSORPTION  LINE 
MEASUREMENTS 

A  computer  program  was  written  to  plot  the  measured 
extinction  coefficient  at  any  CO^  or  DF  laser  frequency  for 
which  FTS  measurements  are  made  as  a  function  of  the  absorp¬ 
tion  coefficient  measured  at  the  peak  of  any  of  several 
atmospheric  H^O  or  HDO  absorption  lines.  The  atmospheric 
line  extinction  coefficient  values  are  taken  to  be  the 
difference  between  the  extinction  coefficient  measured  at 
line  center  and  the  value  measured  some  distance  away  from 
line  center.  It  is  assumed  that  aen  sol  scattering  and 
molecular  continuum  absorption  effects  are  essentially 
constant  over  a  few  cm  ^  and  so  the  difference  between  the 
line  center  and  wing  extinction  coefficient  values  repre¬ 
sents  molecular  abosrption  due  to  that  line.  This  point  is 
discussed  more  thoroughly  in  Reference  2.  Figures  2  and  3 
are  plots  of  extinction  coefficient  at  two  different  DF 
laser  frequencies  versus  molecular  absorption  coefficient 
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FIGURE  2.  EXTINCTION  COEFFICIENT  AT  DF  LASER  LINE  P3(12)  AT  2445.35  cm  IS 
PLOTTED  AS  A  FUNCTION  OF  THE  ABSORPTION  COEFFICIENT  AT  THE  PEAK  OF 
THE  H20  LINE  AT  2844.01  cm-1.  EACH  POINT  REPRESENTS  ONE  DAY.  THE 
CORRELATION  IS  WEAK. 
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at  tnc  peak  cf  the  Ho0  line  at  2844.01  cr.  ~.  Fnrares  4 

L. 

through  8  are  similar  plots  for  5  C0o  laser  frequencies. 
Comparing  the  C09  plots  with  the  DF  plots  is  striking. 

Since  the  water  vapor  peak  absorption  coefficient  is  related 
tc  the  total  water  vapor  content  alone  a  column  in  the  atmos¬ 
phere,  it  appears  that  the  variation  in  extinctance  at  the 
CC>2  frequencies  is  strongly  related  to  variation  in  atmos¬ 
pheric  water  vapor  content,  while  the  variation  at  DF  fre¬ 
quencies  is  due  to  some  other  mechanism.  It  is  still  not 
clear,  however,  whether  the  CO?  dependence  is  due  directly 
to  molecular  absorption  (continuum  effects,  for  instance) 
or  perhaps  to  water  clusters  or  droplets  which  themselves 
depend  on  water  vapor  concentration. 

3.2  SLANT  PATH  TRANSMISSION  MODEL  DEVELOPMENT 

Additional  analysis  of  the  substantial  data  base  of 
slant  path  transmission  spectra  has  proceeded  with  compari¬ 
sons  to  transmission  spectra  predicted  by  our  current  under¬ 
standing  of  atmospheric  absorption.  A  recent  modeling 
approach  developed  by  OptiMetrics  provides  both  a  high 
degree  of  spatial  resolution  of  atmospheric  composition 
with  altitude  and  any  desired  spectral  resolution  up  to 
essentially  infinite  resolution.  A  high  degree  of  spatial 
resolution  is  achieved  by  using  a  25  layer  model  of  the 
earth's  atmosphere  with  variable  layer  spacing.  Between 
the  earth's  surface  and  13  km,  a  minimum  spacing  of  1  km 
is  used;  above  13  km  the  layer  spacing  increases,  reaching 
30  km  for  the  final  layer  between  70  and  100  km.  The 
basic  spectral  resolution  is  infinite  and  ca-  be  convolved 
with  any  desired  "slit  function"  to  provide  comparisons 
with  data  acquired  by  specific  hardware  (as  in  this  case  to 
a  high  resolution  Michelson  interferometer  (FTS)  used  to  ac¬ 
quire  solar  spectra) .  The  spectral  resolution  of  the  model 
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FIGURE  5.  EXTINCTION  COEFFICIENT  AT  C02  LASER  LINE  P(20)  AT  944.19  cm" 1  I 
PLOTTED  AS  A  FUNCTION  OF  THE  ABSORPTION  COEFFICIENT  AT  THE  PEAK 
OF  THE  IlpO  LINE  AT  2844.01  cm"1.  EACH  POINT  REPRESENTS  ONE  DAY 
THE  CORRELATION  IS  RELATIVELY  STRONG. 
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FIGURE  7.  EXTINCTION  COEFFICIENT  AT  THE  C02  LASER  LTNE  R(2R)  AT  J08I.-IK  nil 

IS  PLOTTED  AS  A  FUNCTION  OF  THE  ABSORPTION  COEFFICIENT  AT  THE  FEAR 
OF  THE  H20  LINE  AT  2844.01  cm-1.  EACH  POINT  REPRESENTS  ONF.  DAY  . 
THE  CORRELATION  IS  RELATIVELY  STRONG. 


FIGURE  8.  EXTINCTION  COEFFICIENT  AT  THE  C02  LASER  LINE  R(48)  AT  1091.RR  cm 

IS  PLOTTED  AS  A  FUNCTION  OF  THE  ABSORPTION  COEFFICIENT  AT  THE  PEAK 
OF  THE  110  LINE  AT  2844.(11  cni“  1  .  EACH  pn  |  riT  PEP  RESENTS  nME  DAY. 
THE  CORRELATION  TS  RELATIVELY  STRONG. 
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c  line  by  line  molecular  ansorprior.  calculations, 
description  cf  the  basic  calculations  used  it  the 
v.-iil  be  giver,  here .  The  molecular  absorption  c oof 
k  {  for  the  absorctior.  of  infrared  radiation  is 


=  S  f  (  . 


Yv’he  re 

S  is  the  vibration-rotation  transition  line 
strength  at  temperature  T,  and 

f  (\ )  is  the  normalized  frequency  shape  factor. 

For  near-sea-level  atmospheric  temperature  and  pressure 
the  shape  factor  in  frequency  space  for  each  molecular 
transition  is  closely  approximated  by  the  Lorentz 
function : 


f  (v) 


(l; 


where 

v  is  the  transition  center  frequency,  and 

a  is  the  half  width  at  half  maximum  amplitude 
of  the  absorption  line  profile. 

At  lower  pressure  and  temperature  (i.e.,  higher  altitudes) 
as  collisions  between  molecules  become  more  infrequent, 
the  shape  factor  becomes  a  convolution  of  collision  and 
velocity  broadening  of  each  absorption  transition  process. 
In  these  cases  the  shape  factor  is  computed  as 
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t  15  a r.  integration  parameter. 

This  shape  factor  is  often  referred  to  as  a  Vcict  profile 
and  is  not  normally  used  ir.  line— by-line  caicuiatic.ns  be¬ 
cause  of  its  complexity.  The  Vcict  profile  is  used  in  the 
Opti.Metr ics  slant  path  model  because  of  the  necessity  fer 
high  altitude  calculations.  The  line  strer.oth  £  is  cor¬ 
rected  from  the  291  K  value  given  ir.  the  .-.rc_  line  atlas 
(?_)  according  tc  the  following  approximation : 


where  I  is  the  vibrational  partition  function  which  is 
both  temperature  and  frequency  dependent;  P^T  is  the  value 
of  Pv  at  the  reference  temperature  Tq  (usually  296K) .  The 
rotational  partition  function  PR  is  only  dependent  upon 
temperature;  again  P°  if  the  value  at  the  reference  temp- 
erature.  The  Boltzman  factor  (given  by  the  two  exponentia 
factors  in  Equation  3)  depends  upon  the  lower  state  energy 
Eq  and  the  transition  frequency  vq.  a  more  detailed  treat 
ment  of  this  subject  can  be  found  in  References  5  and  6. 
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S.ar.t  j  it:':  calculations  usin:;  the  Opt i.Mc- tries  model  are 
cased  ei  trier  er.  atmospheric  compos  it  ion  data  contained  ir. 
the  computer  proararr.  for  a  model  atmosphere  or  or.  experi¬ 
mental  input  data  when  available.  The  model  atmospheres 
contain  vertical  density  and  temperature  profiles  fer 
each  absorber  species  which  are  appropriate  for  different 
latitudes  and  seasons  and  are  derive!  from,  an  averaae  c: 
errjoer im.et tal  measurements.  Parameter  values  for  these 
model  at- osjoheres  can  oe  found  in  Reference  5.  Important 
molecular  adsorption  continuum  models  are  included  as 
options  for  nitrogen  and  water  vaoor  in  the  4  to  5  urn  recior. 

Flexibility  cf  comparison  has  been  built  intc  the  slant 
path  transmission  model  through  variable  slant  path  viewing 
angle  and  variable  viewing  platform,  elevations.  The  Voigt 
profile  discussed  earlier  is  incorporated  through  a  hichly 
accurate  and  efficient  computational  subroutine.  Conditions 
warranting  Voigt  profile  calculations  are  tested  internally 
for  each  molecular  transition.  Super-  and  sub-Lorentz  fea¬ 
tures  are  incorporated  where  specific  cases  are  known  to 
occur  as  for  example,  the  sub-Lorentz  far  wings  of  the  4.2 
urn  CO2  bands.  The  transition  line  strengths  and  halfwidths 
are  temperature  corrected  for  accurate  calculations  between 
195  K  and  500  K.  Calculations  outside  this  interval  may 
be  done  at  reduced  accuracy.  The  temperature  dependence 
of  the  1^2  collision  induced  continuum  in  the  3  -  5  urn  recion 
is  that  found  in  LOWTRAN-3B  (7)  and  the  temperature  depen¬ 
dence  of  the  water  vapor  continuum  in  the  8  -  14  ym  region 
is  that  of  Roberts  (8).  Note  should  be  made  here,  that 
calculations  at  wavenumbers  below  800  cm  will  be  slightly 
inaccurate  until  a  rigorous  expansion  of  the  frequency 
dependence  of  the  vibrational  partition  function  is  incor¬ 
porated  into  the  model.  This  can  easily  be  accomplished 
when  necessary  at  the  cost  of  increased  computational  time. 
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2."  ANALYSIS'  OF  MEASURED  FT?  SPECTRA. 

3.3.1  COMPARISONS  OF  SLANT-PATH  TRANSMISSION 
PREDICTIONS  TO  MEASURED  FTS  DATA 

Preliminary  comparisons  of  the  calibrates  slant  path 
"Z2T3.nsrr. x  s  s ion  soo  co  it  c*  *co  c s  1  c u si.  sc.  v  sr^  i ss i "'Ti  cxx,  ^  Yi*-*  s e~ 
viewinc  anaie  were  performed  usina  surface  parameter  cans 
obtained  by  balloon  borne  radiosondes.  The  balloon  radio¬ 
sonde  data  provided  partial  pressure  cf  H-,C,  total  atmos¬ 
pheric  pressure-  and  temperature  up  to  10  km  ir.  altitude. 

Other  molecular  species  concentrations  were  taken  as  mid- 
latitude  summer  model  atmosphere  values  as  were  pressure 
and  temperature  data  when  not  available  from,  balloon  measure¬ 
ments.  Figures  9  through  12  show  comparisons  of  model  pre- 
dictions  (solid  curves)  for  a  39  cm  “  region  near  2500  cm.  ' 

(4  ur,)  .  In  the  cases  considered,  the  White  H00  continuum, 
model  (9)  over-preaicts  absorption  while  the  Burch  1 0 '■ 
model  under-predicted  the  measured  values.  Figure  II  and 

12  show  examples  of  comparisons  of  measured  transmission 

to  calculations  usinc  the  White  H00  continuum  model.  Fiaure 

13  shows  a  comparison  with  the  same  calculation  as  shown  in 
Figure  11  except  that  the  Burch  continuum:  model  is  used. 

It  is  interesting  mo  note  that  the  measured  transmission  is 
unaerpredictec  by  an  amount  in  Figure  II  approximately 
equal  to  the  overprediction  shown  in  Figure  13.  Ficrure  14 
shows  the  agreement  obtained  in  the  10  urn  reuion  usinc  the 
Roberts,  Biberman,  and  Selby  water  continuum;  model  (8). 

All  of  these  calculations  were  done  with  a  20  cm  ^  bound. 
Considering  the  complexity  of  the  25  layer  model  with 
independently  varying  pressure,  temperature  and  molecular 
density  in  each  layer  as  well  as  variable  trace  gas  con¬ 
centrations  and  continuum  absorption  functional  dependence, 
the  agreement  is  remarkable.  It  appears  that  the  continuum 
models  used  here  are  validated  to  a  reasonable  accuracy  and 
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FIGURE  9.  SLANT  PATH  TRANSMISSION  SPECTRA  AT  4  ym  FOR 
87.4°  FROM  ZENITH  WITH  A  1.21  km  SURFACE 
ELEVATION  AND  5.3  TORR  PARTIAL  PRESSURE  H?0 
AT  SURFACE.  THE  SOLID  CURVE  IS  THE  25  LAYER 
MODEL  PREDICTION.  THE  DOTTED  CURVE  IS  THE 
MEASURED  SPECTRA. 
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FIGURE  10.  SLANT  PATH  TRANSMISSION  SPECTRA  AT  4  ym  FOR 
84.4°  FROM  ZENITH  WITH  A  1.23  km  SURFACE 
ELEVATION  AND  5.3  TORR  H20  PARTIAL  PRESSURE 
AT  SURFACE.  THE  SOLID  CURVE  IS  THE  25  LAYER 
MODEL  PREDICTION.  THE  DOTTED  CURVE  IS  THE 
THE  MEASURED  SPECTRA. 
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FIGURE  11.  SLANT  PATH  TRANSMISSION  SPECTRA  AT  4  um  FOR 
78.9°  FROM  ZENITH  WITH  A  1.21  km  SURFACE 
ELEVATION  AND  5.3  TORR  PARTIAL  PRESSURE 
OF  H,0  AT  SURFACE.  THE  SOLID  CURVE  IS  THE 
25  LAYER  MODEL  PREDICTION.  THE  DOTTED  CURVE 
IS  THE  MEASURED  SPECTRA. 
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FIGURE  12. 


SLANT  PATH  TRANSMISSION  SPECTRA  AT  4  Lm  FOR 
78.9°  FROM  ZENITH  WITH  A  1.21  km  SURFACE 
ELEVATION  AND  12.9  TORR  PARTIAL  PRESSURE 
OF  H_0  AT  SURFACE.  THE  SOLID  CURVE  IS  THE 
25  LAYER  MODEL  PREDICTION.  THE  DOTTED  CURVE 
IS  THE  MEASURED  SPICTRA. 
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FIGURE  13.  SLANT  PATH  TRANSMISSION  SPECTRA  AT  4  urn 
WITH  SAME  CONDITIONS  AS  FIGURE  11  EXCEPT 
A  BURCH  WATER  VAPOR  CONTINUUM  IS  USED  INSTEAD 
OF  THE  WHITE  WATER  VAPOR  CONTINUUM. 
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FIGURE  14.  SLANT  PATH  TRANSMISSION  SPECTRA  AT  10.6  yre 

FOR  62.8°  FROM  ZENITH  WITH  A  1.21  km  SURFACE 
ELEVATION  AND  16  TORR  PARTIAL  PRESSURE  OF  H?0 
AT  SURFACE.  THE  SOLID  CURVE  IS  THE  25  LAYER 
MODEL  PREDICTION.  THE  DOTTED  CURVE  IS  THE 
MEASURED  SPECTRA. 
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only  require  minor  adjustments  to  the  temperature  and 
spectral  dependence  of  the  continuum  absorption  coefficient. 
Additional  spectra  should  be  acquired  and  compared  to  the 
code  to  assess  the  temperature  dependence  and  the  varia¬ 
tions  in  spectral  dependence  that  aye  apparent  in  Figure  9. 
Due  to  the  strong  temperature  lapse  in  the  lower  10  km  of 
the  atmosphere  where  significant  amounts  of  occur,  a 

sensitive  test  of  P^O  continuum  absorption  temperature 
dependence  is  now  available.  This  test  is  free  from  problems 
associated  with  instrument  calibrations,  contaminants,  de¬ 
tector  instabilities  and  the  poisoning  or  quenching  effects 
of  limited  laboratory  cell  dimensions  on  long  range  inter- 
molecular  forces  which  may  exist  and  possibly  contribute 
to  continuum  absorption.  The  comparisons  of  the  measured 
FTS  spectra  to  the  slant  path  calculations  for  the  2500  cm  ^ 
to  2530  cm  1  region  shown  in  Figures  9-13  also  demonstrate 

that  the  observed  strengths  of  certain  H ~0  transitions  near 

-1  -1  ^ 

2503  cm  and  2511  cm  disagree  with  the  values  listed  in 

the  AFGL  line  atlas  (4). 

One  of  the  most  important  results  of  this  work  is 
the  capability  to  accurately  predict  slant  path  trans¬ 
mission  at  any  wavelength  contained  in  the  AFGL  line  atlas, 
and  for  any  elevation  angle  and  initial  platform  altitude. 
This  means  that  the  complexity  of  systems'  performance 
analyses  along  atmospheric  slant  paths  can  be  reduced 
to  a  balloon  measurement  of  the  vertical  profile  of  H20 
and  a  computer  calculation.  Of  course  other  spectral  regions 
need  to  be  compared  to  the  AFGL  atlas  and  validated  in  a 
similar  manner  before  confidence  in  the  atlas  values  in 
these  regions  can  be  established.  The  slant  path  trans¬ 
mission  calculation  for  this  25  layer  model  can  be  com¬ 
pleted  in  15  minutes  or  less  depending  upon  the  computer 
used. 
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3.3.2  EXTRACTION  OF  VERTICAL  DENSITY  PROFILES  OF 
ABSORBING  SPECIES  FROM  MEASURED  FTS  DATA 

Another  important  utilization  of  the  calibrated 
slant-path  transmission  spectra  although  much  more  complex 
than  the  comparisons  discussed  thus  far  is  the  extraction 
of  vertical  density  profiles  of  any  IR  absorbing  species 
contained  in  the  line  atlas  from  the  FTS  data.  Using  re¬ 
finements  of  the  comparison  technique  previously  described 
it  becomes  possible  to  substitute  slant  path  solar  FTS 
measurements  for  the  balloon  borne  radiosonde  data.  The 
balloon  monitors  only  those  parameters  for  which  it  is 
specifically  instrumented  —  typically  H20  partial  pressure, 
total  pressure  and  temperature.  Significant  problems  can 
arise  when  the  large  surface  area  of  the  balloon  contami¬ 
nates  the  local  atmosphere  near  the  sensor  upon  ascension 
as  the  balloon  carries  with  it  H20  vapor  from  the  denser 
lower  layers  of  the  atmosphere. 

The  profile  extraction  code  uses  an  algorithm  which 
computes  three  candidate  atmospheric  transmission  spectra 
for  differing  profiles  of  the  species  being  studied  and 
tests  the  root  mean  square  deviation  of  these  test  spectra 
from  the  measured  solar  spectra  being  fit.  A  best  profile 
is  selected  and  the  procedure  is  repeated  using  the  new 
best  atmospheric  concentration  profile  as  a  mean  and 
smaller  deviations  than  in  the  previous  iteration.  Im¬ 
proved  convergence  to  the  best  rms  fit  was  achieved  by 
scaling  the  first  10  km  path  segment  without  a  vertical 
gradient  at  first  and  then  fitting  a  skew  factor  to  this 
10  km  segment  before  stepping  through  each  layer  indivi¬ 
dually,  starting  at  the  surface,  to  look  for  small  irregu¬ 
larities.  This  approach  was  found  to  be  satisfactory  for 
all  cases  except  possible  large  temperature  inversions, 
in  which  case  an  average  result  was  obtained  that  did  not 
reproduce  the  inversion. 
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Figures  15a  through  19a,  show  the  results  of  the  best 
fit  computed  spectrum  compared  to  the  measured  solar  spectra 
and  Figures  15b  through  19b  show  the  corresponding  vertical 
H2O  concentration  profiles  compared  to  balloon  data  taken 
close  in  time.  The  comparisons  show  generally  good  agree¬ 
ment  near  the  ground  with  the  balloon  reading  showing 
relatively  higher  water  vapor  at  middle  elevations  and 
the  two  profiles  in  much  better  agreement  above  7  km.  It 
should  be  pointed  out  again  that  a  possible  explanation  for 
the  balloon  reading  being  high  on  ascension  even  when  cali¬ 
brated  to  an  absolute  reading  at  the  surface  is  that  water 
vapor  absorbed  on  the  large  surface  area  of  the  balloon 
is  carried  from  regions  of  higher  to  lower  density  as  the 
balloon  climbs.  If  we  accept  the  balloon  data  as  accurate, 
then  the  model  values  for  the  continuum  should  be  ad¬ 

justed  to  produce  better  agreement. 

The  uncertainty  in  the  extracted  vertical  density 
profiles  was  checked  by  choosing  a-priori  an  arbitrary 
water  density  profile  and  inputting  this  profile  into  the 
25  layer  transmission  model  to  generate  a  synthetic  slant- 
path  transmission.  This  noise-free  spectrum  was  then  used 
by  the  25  layer  vertical  density  profile  fitting  routine 
to  extract  a  profile.  The  results  of  this  test  case  are 
shown  in  Figure  20a.  Here  the  solid  curve  is  the  trans¬ 
mission  spectrum  corresponding  to  the  assumed  ^0  density 
profile  and  the  dotted  curve  is  the  least ' squares  fit  to 
this  test  spectrum  obtained  with  25  layer  model  using  the 
vertical  density  profile  of  H2O  as  the  parameter  being 
fitted.  Figure  20b  shows  the  H2O  density  profile  resulting 
from  the  fit  to  the  synthetic  transmission  spectrum  with  no 
other  inputs.  As  can  be  seen,  maximum  deviations  in  H2O 
density  with  altitude  are  less  than  10%  with  most  points 
much  closer. 
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FIGURE  15a.  THE  SOLID  CURVE  IS  A  LEAST  SQUARES  FIT  OF  THE 

25  LAYER  MODEL,  WITH  H20  AS  THE  SEARCH  PARAMETER, 
TO  THE  MEASURED  SLANT  PATH  TRANSMISSION  SPECTRUM 
SHOWN  AS  THE  DOTTED  CURVE  FOR  THE  CONDITIONS 
GIVEN  ON  THE  FIGURE. 
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FIGURE  15b.  THE  H20  VERTICAL  DENSITY  PROFILE  RESULTING  FROM 
THE  BEST  FIT  SPECTRUM  GIVEN  IN  FIGURE  15a  IS 
SHOWN  AS  THE  DOTTED  CURVE.  THE  SOLID  CURVE  IS 
A  BALLOON  BORNE  RADIOSONDE  MEASUREMENT  TAKEN 
THE  SAME  DAY  IN  THE  MORNING. 
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FIGURE  16a.  THE  SOLID  CURVE  IS  A  LEAST  SQUARES  FIT  OF  THE 
25  LAYER  MODEL,  WITH  H?0  AS  THE  SEARCH  PARA¬ 
METER,  TO  THE  MEASURED  SLANT  PATH  TRANSMISSION 
SPECTRUM  SHOWN  AS  THE  DOTTED  CURVE  FOP  THE 
CONDITIONS  ON  THE  FIGURE. 
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FIGURE  16b.  THE  H20  VERTICAL  DENSITY  PROFILE  RESULTING 
FROM  THE  BEST  FIT  SPECTRUM  GIVEN  IN  FIGURE 
16a  IS  SHOWN  AS  THE  DOTTED  CURVE.  THE  SOLID 
CURVE  IS  A  BALLOON  BORNE  RADIOSONDE  MEASURE¬ 
MENT  TAKEN  THE  SAME  DAY  IN  THE  MORNING. 
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FIGURE  17a.  THE  SOLID  CURVE  IS  A  LEAST  SQUARES  FIT  OF  THE 
25  LAYER  MODEL,  WITH  H-0  AS  THE  SEARCH  PARA¬ 
METER,  TO  THE  MEASURED  SLANT  PATH  TRANSMISSION 
SPECTRUM  CONDITIONS  GIVEN  ON  THE  FIGURE. 


36 


ALTITUDE  (KM) 


H20  VERTICAL  PROFILE 


SURF  ELEV.  <KM>  1. 21 
SURF  TEMP.  <F>  81.5 

SURF  BAR  PR  <TR>  659.  1 
LOOK  ANGLE  <DEG>  62.  B 


_  BALLOON  DATA 

TESTAE  SOLAR  SPECT. 
.  .  .  FTS  FIT  TO  SOLAR 


FIGURE  17b.  THE  H20  VERTICAL  DENSITY  PROFILE  RESULTING 

FROM  THE  BEST  FIT  SPECTRUM  GIVEN  IN  FIGURE  17a 
IS  SHOWN  AS  THE  DOTTED  CURVE.  THE  SOLID  CURVE 
IS  A  BALLOON  BORNE  RADIOSONDE  MEASUREMENT 
TAKEN  THE  SAME  DAY  IN  THE  MORNING. 
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FIGURE  18a.  THE  SOLID  CURVE  IS  A  LEAST  SQUARES  FIT  OF  THE 
25  LAYER  MODEL,  WITH  H20  AS  THE  SEARCH  PARA¬ 
METER,  TO  THE  MEASURED  SLANT  PATH  TRANSMISSION 
SPECTRUM  SHOWN  AS  THE  DOTTED  CURVE  FOR  THE 
CONDITIONS  GIVEN  ON  THE  FIGURE. 
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FIGURE  18b.  THE  H20  VERTICAL  DENSITY  PROFILE  RESULTING 

FROM  THE  BEST  FIT  SPECTRUM  GIVEN  IN  FIGURE  18a 
IS  SHOWN  AS  THE  DOTTED  CURVE.  THE  SOLID  CURVE 
IS  A  BALLOON  BORNE  RADIOSONDE  MEASUREMENT  TAKEN 
THE  SAME  DAY  IN  THE  MORNING. 
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FIGURE  19a. 


THE  SOLID  CURVE  IS  A  LEAST  SQUARES  FIT  OF 
THE  25  LAYER  MODEL,  WITH  H,0  AS  THE  SEARCH 
PARAMETER,  TO  THE  MEASURED  SLANT  PATH 
TRANSMISSION  SPECTRUM  SHOWN  AS  THE  DOTTED 
CURVE  FOR  THE  CONDITIONS  GIVEN  ON  THE  FIGURE. 
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FIGURE  19b. 


THE  H-0  VERTICAL  DENSITY  PROFILE  RESULTING 
FROM  THE  BEST  FIT  SPECTRUM  GIVEN  IN  FIGURE 
18a  IS  SHOWN  AS  THE  DOTTED  CURVE.  THE 
SOLID  CURVE  IS  A  BALLOON  BORNE  RADIOSONDE 
MEASUREMENT  TAKEN  THE  SAME  DAY  IN  THE 
MORNING. 


41 


TRANSMISSION 


SLANT  PATH  TRANSMISSION 


SURFACE  LAYER.* 

H20  15.  00  CTPRR5 

C02  .  212 

03  .  224E-04 

N20  .  1B0E-03 

CO  .  48 1 E  — 04 

CH4  .  103E-E2 

02  134.  4 

LOOK  ANGLE  (DEO  90.  0 
SURF  ELEV.  (KM)  1. 2 
SURF  TEMP  (F>  57.  8 

SURF  BAR  PR  (TR) 64 1.4 
MIDLAT  SUMMER  ATM. 
BURCH  N2  CONTINUUM 
WHITE  H2D  CONTINUUM 
BOUND  (CM- 15  5. 

_  SYNTHETIC  PROFILE 

.  ..  25  LAYER  FIT  SPECT. 


FIGURE  20a.  THE  SOLID  CURVE  IS  THE  SLANT  PATH  TRANSMISSION 
GENERATED  BY  A  CHOSEN  H20  DENSITY  VERSUS 
ALTITUDE  PROFILE.  THE  DOTTED  CURVE  IS  THE 
RESULT  OF  A  25  LAYER  LEAST  SQUARES  FIT  TO 
THE  TRANSMISSION  PROFILE  (WITH  NO  KNOWLEDGE 
OF  ASSUMED  H20  PROFILE) . 
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FIGURE  20b. 


THE  SOLID  CURVE  IS  THE  ASSUMED  H20  DENSITY 
VERSUS  ALTITUDE  PROFILE.  THE  DOTTED  CURVE 
IS  THE  H20  DENSITY  VERSUS  ALTITUDE  EXTRACTED 
BY  THE  25  LAYER  FIT  ROUTINE  BASED  SOLELY  ON 
THE  TRANSMISSION  SPECTRUM. 


This  test  case  is  a  very  useful  method  for  checking 
the  validity  of  the  procedure  outlined  above  for  extracting 
the  vertical  density  profile  of  a  particular  molecular 
absorber.  Since  no  noise  is  present  in  the  "ideal"  syn¬ 
thetic  spectrum  which  contains  a  "known"  vertical  density 
distribution  of  the  molecular  absorber  under  study,  and 
since  no  assumptions  regarding  the  accuracy  of  spectral 
line  parameters  and  continuum  absorption  coefficients  are 
involved,  the  comparisons  shown  in  Figure  20b  tests  the 
computational  procedure  alone,  independent  of  the  addi¬ 
tional  factors  involed  in  extraction  of  profiles  from  the 
FTS  data  and  comparisons  to  the  balloon  sonde  data. 

An  improved  fitting  routine  has  been  set  up  but  not 
tested  which  will  use  the  unique  temperature  of  each  layer 
as  a  discriminant  of  the  distribution  of  the  total  burden 
of  the  molecule  under  study  should  be  distributed.  The 
uniqueness  of  the  differing  ground  state  energy  levels 
and  thermal  populations  for  the  various  transitions  is 
included  in  the  molecular  density  distribution  fitting 
procedure . 

The  calibrated  solar  spectra  collected  during  the 
past  two  years  provide  a  statistically  significant  data 
base  of  surface  concentrations,  vertical  density  profiles, 
and  variations  in  the  path  integrated  effects  of  the  primary 
IR  absorbing  molecules  and  their  isotopes.  The  most  im¬ 
portant  one  of  these  in  the  chemical  lasci  region  is  HDO 
which  is  difficult  to  measure  directly  by  means  other 
than  the  Fourier  Transform  Spectroscopy  used  here. 

3.3.3  FURTHER  UTILIZATION  OF  EXISTING  SLANT  PATH 
FTS  SOLAR  SPECTRA 

The  problem  of  modeling  long  range  detection  of  infra¬ 
red  target  signatures  can  be  addressed  by  use  of  FTS  solar 
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spectra  of  low  elevation  angle  paths.  Transmission  spectra 
at  angles  near  1°  above  the  horizon  have  been  measured  and 
provide  useful  data  for  long  path  detection  analyses.  Figure 
9  shows  an  example  of  long  path  transmission  and  fair  agree¬ 
ment  with  the  model.  The  model  computation  for  this  case 
can  be  evaluated  for  the  relative  importance  of  the  various 
absorption  components  at  any  specified  frequency  providing 
a  better  understanding  of  the  problem. 

Another  application  involving  potential  utilization  of 
the  slant  path  solar  data  is  the  detection  and  identifica¬ 
tion  of  trace  and  pollutant  gases  which  absorb  in  the  in¬ 
frared.  Information  concerning  abnormally  high  or  low 
molecular  concentrations  and  concentration  fluctuations  of 
the  more  abundant  atmospheric  molecules  is  also  an  important 
problem  which  can  be  addressed  using  the  solar  spectra. 
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4 

CONCLUSIONS  AND  RECOMMENDATIONS 

Upgraded  FTS  operating  system  software  has  been  im¬ 
plemented  on  the  system  being  used  for  slant  path  solar 
spectral  measuements  and  tested  as  discussed  in  Sections 
2.1  and  2.3.  The  new  software  allows  increased  spectral 
resolution  to  be  used,  facilitating  critical  comparisons 
of  measured  spectra  to  calculations  of  slant  path  atmos¬ 
pheric  transmission.  Certain  modifications  in  data  collec¬ 
tion  procedures  were  required  to  handle  the  larger  data 
sets,  but  these  have  been  satisfactorily  accomplished  so 
that  the  former  data  collection  rate  can  be  maintained. 

The  FTS  system  installation  in  the  field  measurement 
van  has  been  completed  allowing  continued  collection  of 
slant  path  solar  spectra  to  be  carried  out  from  a  variety 
of  measurement  sites.  In  this  way  local  influences  peculiar 
to  any  of  these  sites  may  be  studied.  Although  the  major 
segment  of  the  atmospheric  slant  path  over  which  the  trans¬ 
mission  measurement  is  made  is  essentially  independent  of 
any  particular  location  of  the  FTS  system,  the  dominant 
influence  of  the  atmospheric  surface  layer  (^10  km)  on 
absorbing  gas  composition  and  density  allows  "site  peculiar" 
contributions  to  total  absorption  along  the  slant  path  to  be 
studied. 

Extinction  coefficients  at  selected  DF  and  CO^  laser 
lines  were  derived  from  several  examples  of  the  FTS  spectra 
and  plotted  against  measured  peak  water  vapor  absorption 
coefficients  for  the  1^0  line  at  2844.01  cm  ^  obtained 
from  the  same  spectra.  The  results  of  these  comparisons 
are  shown  in  Figures  2-8.  The  correlation  of  the  DF 
P^(12)  and  P^(7)  line  data  (Figures  2  and  3)  with  measured 
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H^O  absorption  coefficients  is  weak  and  relatively  independ¬ 
ent  of  increasing  absorption,  indicating  the  relative 

importance  of  extinction  mechanisms  other  than  water  vapor 
at  the  frequencies  of  the  two  DF  laser  lines.  Both  of 
these  lines  are  located  in  an  ^0  absorption  band.  In 
addition  the  N2  pressure-induced  continuum  absorption  also 
contributes  at  these  frequencies.  By  contrast  the  C02 
laser  line  data  shown  in  Figures  4-8  show  good  correla¬ 
tion  with  measured  H20  absorption  coefficients  demonstrat¬ 
ing  that  water  vapor  absorption  is  the  dominant  extinction 
mechanism  at  these  frequencies. 

The  slant  path  transmission  model  development  described 
in  Section  3.2  represents  significant  progress  in  the  utili¬ 
zation  of  the  slant  path  solar  spectra  collected  during  the 
past  2-1/2  years.  Variability  of  a  large  number  of  modeling 
parameters  may  now  be  studied  by  comparing  previously  col¬ 
lected  FTS  data  with  the  slant  path  model  predictions. 
Improved  understanding  of  ambient  atmospheric  gas  composi¬ 
tions  and  their  variability  as  well  as  identification  of 
possible  discrepancies  in  the  AFGL  line  parameter  listing 
should  now  be  possible.  Results  of  initial  FTS-slant  path 
modeling  comparisons  are  shown  in  Figures  9-13  (Section 
3.3.1).  Initial  indications  based  on  these  comparisons  show 
that  the  3  -  5  ym  water  vapor  continuum  absorption  falls  be¬ 
tween  the  two  competing  models,  at  least  for  the  atmospheric 
conditions  encountered  and  for  the  spectral  interval  exa¬ 
mined  (2500  cm-1  to  2530  cm-'*")  .  The  above  comparisons  also 
point  out  disagreements  between  observed  line  strengths  and 
computed  values  using  the  most  current  AFGL  parameters  for 
water  vapor  absorption  lines  near  2503  cm  ^  and  2511  cm 
The  comparison  of  FTS  measurements  to  slant  path  calcula¬ 
tions  in  the  CC>2  laser  region  between  525  cm  "*"  and  955  cm  ^ 
that  is  shown  in  Figure  14  indicates  excellent  agreement. 
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No  modeling  discrepancies  are  apparent  in  the  data  shown 
in  Figure  14,  however  additional  comparisons  must  be 
carried  out  before  conclusive  results  can  be  obtained. 

Use  of  the  slant  path  modeling  technique  to  extract 
vertical  profiles  of  particular  molecular  absorbers  is 
described  in  3.3.2.  The  dominant  atmospheric  absorber, 
in  the  spectral  region  studied,  namely  water  vapor  is  used 
as  an  example.  The  selection  of  test  criteria  and  hier¬ 
archy  of  steps  involved  in  the  profile  extraction  procedures 
are  illustrated  by  means  of  several  examples  shown  in  Figures 
15  -  19.  In  these  figures  comparisons  between  vertical  J^O 
density  profiles  extracted  from  the  FTS  data  to  balloon 
sonde  data  are  shown.  Good  results  were  obtained  in  the 
comparisons  however  certain  limitations  arise  in  the  appli¬ 
cation  of  these  techniques  when  pronounced  temperature  (and 
consequently  water  vapor  density)  inversions  occur  as  a 
function  of  altitude. 

A  test  of  the  uncertainty  associated  with  the  profile 
extraction  procedure  is  described  and  results  obtained  in 
test  calculations  are  presented  in  F’igure  20.  The  test 
showed  the  ability  of  the  extraction  procedure  to  faith¬ 
fully  reproduce  an  I^O  profile  which  was  used  to  calculate 
a  synthetic  test  spectrum  which  was  then  substituted  for  an 
experimental  slant  path  spectrum  in  the  extraction  proce¬ 
dure.  Further  development  of  this  technique  will  be  pur- 
sured  since  it  establishes  a  procedure  for  using  the  FTS 
data  in  a  "stand-alone"  mode  without  the  requirement  for 
simultaneous  atmospheric  soundings  for  temperature,  water 
vapor  density,  etc.  As  discussed  in  Section  3.3.2,  uncer¬ 
tainties  caused  by  atmospheric  contamination  from  the  bal¬ 
loon  platform  itself  are  avoided  when  the  FTS  spectra  are 
used  alone. 


48 


REFERENCES 


1.  B.L.  Bean,  W.L.  Flowers,  W.M.  Gutman,  A.V.  Jelinek 
and  R.L.  Spellicy,  Laser,  IR  and  NMMW  Propagation 
Measurements  and  Analysers'/  OptiMetrics ,  Inc.  Report 
OMI- 80-001 ,  1980. 

2.  W.  Flowers,  W.  Gutman,  Atmospheric  Characterization 
Measurements  at  White  Sands  Missile  Range;  Point  ~ 
Sampling  and  Integrated  Path  Results,  SAI  Report 
SAI-79-002-AA,  1979. 

3.  W.  Gutman,  HEL  Related  Gas  Analysis:  Part  II:  Fourier 
Transform  Solar  Spectrum  Measurements,  SAI  Report 
SAI-79-007-AA-II ,  1979. 

4.  R. A.  McClatchey,  et  al.,  AFCRL  Atmospheric  Absorption 
Line  Parameters  Compilation,  Air  Force  Research 
Laboratories,  L.H.  Hanscom  Field,  Bedford,  MA, 

AFCRL- TR- 73-0096,  1973;  L.S.  Rothman,  "Update  of  the 
AFGL  Atmospheric  Absorption  Line  Parameters  Compila¬ 
tion1/  Appl .  Opt.  17 ,  3517,  1978. 

5.  W.L.  Wolfe  and  G.J.  Zissis,  The  Infrared  Handbook, 
IRIA,  Environmental  Research  Institute  of  Michigan, 
1978,  Chapter  5. 

6.  F.S.  Mills,  Absorption  of  Deuterium  Fluoride  Laser 
Radiation  by  the  Atmosphere,  Ohio  State  Electro- 
Science  Laboratory  Report  4054-3,  1975. 

7.  J.E.A.  Selby,  E.P.  Shettle  and  R. A.  McClatchey, 
Atmospheric  Transmittance  from  .25  to  .28  ym; 
Supplement  LOWTRAN  3B  (1976) ,  Air  Force  Geophysics 
Laboratory,  Bedford,  MA,  AFGL-TR-76-0258 ,  1976. 

8.  R.E.  Roberts,  L.M.  Biberman  and  J.E.A.  Selby, 

Infrared  Continuum  Absorption  of  Atmospheric  Water 
Vapor  m  the  8-12  ym  Window,  Paper  P-1184,  Insti¬ 
tute  for  Defense  Analyses,  Arlington,  VA,  1976. 

9.  W.R.  Watkins  and  K.O.  White,  Water  Vapor  Continuum 

Absorption  Measurements  (3,5-  4.0  ym)  Using  HDO 
Depleted  Water,  Optics  Letters,  1^,  31  ,  (1977). 

10.  D.E.  Burch,  D.A.  Grynak,  and  J.D.  Pembrook,  Investi¬ 
gation  of  the  Infrared  Radiation  by  Atmospheric  Gases; 
Water,  Nitrogen,  Nitrous  Oxide,  Philco-Ford  Aeronu- 
tronic  Division  Report  U-4784,  AFCRL-71-0124 ,  1974. 


49 


